Quantum dynamics of vortices in mesoscopic magnetic disks 
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Model of quantum depinning of magnetic vortex cores from line defects in a disk geometry and 
under the application of an in-plane magnetic field has been developed within the framework of the 
Caldeira-Leggett theory. The corresponding instanton solutions are computed for several values of 
the magnetic field. Expressions for the crossover temperature Tc and for the depinning rate r(r) 
are obtained. Fitting of the theory parameters to experimental data is also presented. 
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I. INTRODUCTION 

Macroscopic quantum tunneling of mesoscopic solid- 
state objects has been intensively studied in the past. Ex- 
amples include single domain particlesi^— , domain walls 
in magnets^"—, magnetic clustersii^, flux lines in type-II 
superconductors^iiS and normal-superconducting inter- 
faces in type-I superconductor a ^ . It is well known 
that micron-size circular disks made of soft ferromag- 
netic materials exhibit the vortex state as the ground 
state of the system for a wide variety of diameters and 
thicknesses^^—. This essentially non-uniform magnetic 
configuration is characterized by the curling of the mag- 
netization in the plane of the disk, leaving virtually 
no magnetic "charges" . The very weak uncompensated 
magnetic moment of the disk sticks out of a small area 
confined to the vortex core (VC). The diameter of the 
core is comparable to the material exchange lengt h^^d^ 
and, because of the strong exchange interaction among 
the out-of-plane spins in the VC, it behaves as an inde- 
pendent entity of mesoscopic size. 

Recent experimental works have reported that the dy- 
namics of the VC can be affected by the presence of struc- 
tural defects in the samplei^"— . This is indicative of the 
elastic nature of the VC line, whose finite elasticity is 
provided by the exchange interactionSS. In Ref. [l^ non- 
thermal magnetic relaxations under the application of an 
in-plane magnetic field are reported below T = 9 K. It is 
attributed to the macroscopic quantum tunneling of the 
elastic VC line through pinning barriers when relaxing 
towards its equilibrium position. In such range of low 
temperatures only the softest dynamical mode can be 
activated, which corresponds to the gyrotropic motion of 
the vortex state. It consists of the spiral-like precessional 
motion of the VC as a whole^^"— and it is intrinsically 
distinct from conventional spin wave excitations. It can 
also be viewed as the uniform precession of the magnetic 
moment of the disk due to the vortex. 

The aim of this paper is to study the mechanism of 
quantum tunneling of the elastic VC line through a pin- 
ning barrier during the gyrotropic motion. We focus our 
attention on line defects, which can be originated for in- 



stance by linear dislocations along the disk symmetry 
axis. This case may be relevant to experiments performed 
in Ref. [2^ since linear defects provide the maximum pin- 
ning and, therefore, the VC line in the equilibrium state 
is likely to align locally with these defects. Such a sit- 
uation would be similar to pinning of domain walls by 
interfaces and grain boundaries. Thus, we are consid- 
ering the depinning of a small segment of the VC line 
from a line defect. The problem of quantum and thermal 
depinning of a massive elastic string trapped in a linear 
defect and subject to a small driving force was considered 
by Skvortso\*2S. The problem studied here is different as 
it involves gyrotropic motion of a massless vortex that is 
equivalent to the motion of a trapped charged string in 
a magnetic field^^. We study this problem with account 
of Caldeira-Leggett type dissipation. 

The paper is structured as follows. In Sec. II the 
Lagrangian formalism of the generalized Thiele's equa- 
tion is presented and Caldeira-Leggett theory is applied 
to obtain the depinning rate. The imaginary-time dy- 
namical equation for instantons is derived in Section III 
and numerical solutions are computed. In Section IV the 
crossover temperature between the quantum and thermal 
regime is obtained. Discussion and fitting of the theory 
parameters (which is related to the pinning potential) 
to experimental data are provided in Sec. V. Also final 
conclusions are included in this section. 



II. ELASTIC THIELE'S LAGRANGIAN 
FORMALISM AND DEPINNING RATE 

In this paper we restrict ourselves to a circular disk 
geometry and to an applied in-plane magnetic field con- 
figuration. The VC line is pinned by the line defect going 
in the Z direction (symmetry axis of the disk) at the cen- 
ter of the disk. The vortex line shall be described by the 
vector field X ~ {x,y), where x{t,z) and y{t,z) are co- 
ordinates of the center of the VC in the XY plane. The 
dependence on the Z-coordinate emerges from the elas- 
tic nature of this magnetic structure. Figure [1] shows an 
sketch of the vortex line deformation due to pinning and 
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its gyroscopic motion. 




FIG. 1: (Color online) Vortex state and depinning via nucle- 
ation of the part of the VC line in a circular disk made of soft 
ferromagnetic material. 

The softest dynamical mode of the VC, and hence of 
the whole vortex, originates from gyroscopic motion and 
it is described by the generalized Thiele's equation^^: 



X{t,z) X pG + dJl^ + V^uj = 



(1) 



time derivative. The gyrovector 
pce-z is responsible for the gyroscopic 



where ' ' ' means 
density2£ pQ 

motion of the VC and its modulus is given by pQ — 
27rprti,Afs/7, where Ms is the saturation magnetization, 
7 is the gyromagnetic ratio, p — ±1 is the polariza- 
tion of the VC and riy — ±1 is the vorticity of the 
magnetization of the disk. The potential energy den- 
sity u){X,dzX) spUts into the sum of two contributions, 
ijJi{X) and ijJ2{dzX). The latter is the elastic energy 



term, Lxj2{dzX) = 



— — I , which is provided by the 

exchange interaction. The elastic constant is given by 
A = 27rA ln(i?/Ao), where R is the radius of the disk, 
A is the exchange constant and Ao = \J AjM^ is the 
exchange length of the ferromagnetic material. Finally, 
= —5ijj / 5{dzX) = —XdzX is the generalized momen- 
tum density with respect to Z. Consequently, the gener- 
alized Thiele's equation becomes 



X{t, z)x pG- \dlX{t, z) + V^w = 



(2) 



Let L be the thickness of the circular disk. The La- 
grangian corresponding to the above equation is given 



c[t,x,x,dzX] 



Az IX ■ An 



uj(X,dzX) 



(3) 



where Ap^ = pcysx is the gyrovector potential in a con- 
venient gauge^fl. The depinning rate at a temperature T, 
r(r) = A{T) exp [-B{T)], is obtained by performing the 
imaginary-time path integral^ 



J D{x} J i^{y}exp -1 AtCe 



(4) 



over Xr = X{t, z) trajectories, which are periodic in r 
with period h/ksT. Notice that t — it is the imaginary 
time and Ce is the Euclidean version of Eq. ([3]). That 
is. 



CE[T,XT,XrjdzXr] = J dz iX-r ■ Ap^ 

+ ijiXr,dzXr) 



(5) 



The energy density u;i{Xr) splits into the sum of three 
terms: The first one, uxviXr), represents the sum of the 
magnetostatic and exchange contributions in the z-cross- 
section, whose dependence on the vortex core coordinates 
is LUxviXr) ^ X^ for small displacements^^. The second 
term, ajdepiX.,-), represents the pinning energy density 
associated to the line defect. Recent experimental works 
have reported an even quartic dependence of pinning po- 
tentials on the VC coordinates for small displacements in 
permalloy rings^S. Consequently, it is legitimate to take 
the following functional dependence for the sum of both 
terms: 

{UJXY + (^dep)iXr) = {x^ + y^) - ^ /3p4 (x^ , ) (6) 

where (k, /?) are the parameters of our model and Pi^x, y) 
is a linear combination of monomials of degree four on 
variables x and y. The last term is the Zeeman energy 



z X Hi, 



■Xr 



density, which is given by^ ujziXr) = —p. 

-with n = (2TT/3)MsnyR- for small displacements. The 
latter correspond to the application of a weak in-plane 
magnetic field Hin- In what follows, Hin = ^Hcy is 
applied along the Y direction. 

The simple dependence pi{xr, J/t) = keeps the main 
features of the pinning potential (see Section V). We also 

neglect the elastic term ^A {^§f^ ■ From all these con- 
siderations, the Lagrangian ([SJ becomes 



CE[T,Xr,Xr,dzXr]= I Azl -ipayrXr-phXr 



) (7) 



2~r . 4"^^ ■ 2 V dz 

Finally, Gaussian integration over y^. reduces Eq. (|3|) 
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to 



with 



D{x} exp 



drC 



E,eff 



(8) 



2 \ K 



K 2 Pi ^ f dXr 



(9) 



Within the framework of the Caldeira-Leggett 
theory^, dissipation is taken into account by adding a 
term 

^ [\z{dr[ ^,^ i-rir,z)-xAn,z)? 

to the action of Eq. (|8]). The dissipative constant -q is 
related to the damping of the magnetic vortex core^ and 
Ref. m shows that rj ~ SaLLclPcl, with Ullg being the 
Gilbert damping parameter. Introducing dimensionless 
variables f = {k/^/2\pg\)t, z = [nj^X)^!'^ z and u = 
(2j3 / k^/'^Xt, the depinning exponent becomes 



B{T, h) 



2h/3 

V 

Pg\ 



2V27r| 



dz d> dr 



dfi 



{u{t,z) - u{Ti,z)y 

(f-fl)2 



(11) 



where ' ' ' means derivative with respect to z, V{u, h) = 
— hu + u 4~ ^^'^ normalized energy potential and 

h — 2^2(5 / fiH . Let uo{h) be the relative minimum of 
V for a fixed value of h. We reescale the energy potential 
V{u, h) V{u, h) := V{uoih) + u,h)- V{uo{h), h) and 
the variable u — ^ uo(/i) + it so that we obtain V{u, h) — 



"oC*) ) ~ uo{h)u — -u' 



III. INSTANTONS OF THE DISSIPATIVE 1+1 
MODEL 

Quantum depinning of the VC line is given by the in- 
stanton solution of the Euler-Lagrange equations of mo- 
tion of the 1+1 field theory described by Eq. ([TT|). This 
gives 

ii + -u" - (2 - 'iul{h))u + 3uo{h)u'^ + u^- 

/dn "("';^-"(J'"^ =0 (12) 

TT \pg\ Jr (t-ti)2 
with boundary conditions 

u{-n/2,z) ^u{n/2,z) zeR 

max u(f, z) = u(0, z) zgR (13) 
Te[-n/2,n/2] 



that must be periodic on the imaginary time r with the 

period £7 = — ^= . This equation cannot be solved 

V2\pg\ ksT 

analytically, so we must proceed by means of numerical 
methods. Notice that in the computation of instantons 
we can safely extend the integration over z in Eq. (jlip 
on the the whole set of real numbers. 



A. Zero temperature 

In this case we apply the 2D Fourier transform 

u{uj,9)^^ [ dfdzu(f,f)e*("^+''^") (14) 
2tt ./b2 



to Eq. ((T2|) and obtain 

u{uj,e) 



1 / 3uo(fe) 

Lj^ + 0^ + V2\^\T]/\pG\ + 2~3ul{h)\ 2tt 

1 

(2^)2 



X / dLL!id9iu{uji,0i)u{uj — u}i,0 — 9i) + -j^—^x 



d^uj d^0 u{liJ2, 92)u{u!i—U!2, 0l—02)u{uj~iOi, 0~0i) 



(15) 

which is an integral equation for u. The depinning expo- 
nent (fTTj) in the Fourier space becomes 



B{T = 0,/i) = 
3 



\pg\Vk\ 



du d0 u{uj,0)u{-uj,-0) 



2h(3 

^2 + 6*2 |a;| 7] \ uo{h) 



u{uJi,0i)u{uj2,02)H-UJl-UJ2, -01-02) 



2-K 
1 

(4^ 



d^wd^ 



u{l^l,0l)u{l^2,02)u{l^3,0:i)u{-UJl-UJ2-iOj., -01-02-0?.) 



(16) 



The zero-temperature instanton is computed using an 
algorithm that is a field-theory extension of the algorithm 
introduced in Refs. [3^ [35! for the problem of dissipative 
quantum tunneling of a particle: To begin with, we in- 
troduce the operator 



0(A, a, u(w, 0),h) = 



1 



w2 + 02 + V2|c^|?//Ipg| + 2 - iul{h) 
\ I du)id0iu{uji,0i)u{ijj — uJi,0 — 0i)+ 



a I d cD d u{uj2,02)u{uJi—UJ2,0l—02)u{uj—LjJl,0—0l) 



(17) 
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Secondly, it is important to point out the scaling prop- 
erty of this operator because it will be used in the com- 
putation of Eq. (fT6| : Given any triplet (Aq, ao, Mo I ^)) 
satisfying the identity P3|) . so will any other triplet 
{Xi,ai,ui{ui,0)) provided that 



"1 = ao/x^> 



(18) 
(19) 
(20) 



where x is a constant. This means that if we are able to 
find a solution (Ai, ai, ui(w, 9)) for arbitrary parameters 
(Ai,ai), then we can obtain the solution corresponding 
to the pair (Aq, ao) simply by rescaling ui{uj, 9) by a fac- 
tor X = Ai/Ao as long as (Ai/Ao)^ = ai/ao is verified. 
The algorithm consists of the following steps: 

1. Start with an initial (Aq, ao, 'Uo('^: ^))- 

2. Let Ui{uj, 9) — 0(Ao, ao, uq{u!, 9), h). 

3. Calculate Ai = Xq/x^,C(i — o-o/x^j where x = 
mi(0)/mo(0). 

4. Find ^2(^,6') = 0{Xi,ai,ui{u},9),h). 

5. Repeat steps (2)-(4) until the successive difference 
satisfies a preset convergence criterion. 

The output is the triplet (A„, a„, u„(w, 6*)). The final 
step consists of reescaling u„ to obtain the solution corre- 
sponding to the pair (A, a) — {3uo{h)/27r, l/(27r)^): from 
the scaling property we know that the reescaling rules of 
the A- and a- terms of Eq. (|15|) arc different. Thus, to 
obtain an accurate approximation of the instanton solu- 
tion we have split u{uj^ 9) into the sum of two functions 
ui{ijj,9) and U2{ijJ,9) in the above algorithm, and calcu- 
lated their next iteration by means of the A-term, respec- 
tively the a-term of the operator (IT7)) . Finally, we rescale 
Ml by a factor 27rA„/3wo(/i) and 112 by a factor 
The depinning rate is calculated evaluating Eq. (|T5l) at 
this solution. 



B. Non-zero temperature 

In the r 7^ case, taking into account the finite peri- 
odicity on f we consider a solution of the type 



x(f,z) = ^u„(z)e-*""^ 



(21) 



27r7i 

with w„ = ^ for all n e Z. Introducing this functional 
dependence into Eq. (jl2l) and applying a ID Fourier 



transform we obtain 

Un{9) 



-y 

27r ^ 



ujI + 9^ + V2|w„|77/Ipg| + 2 - iulih) 
(iu^y^f ^e,u,{9,)u^_,{9-9,)+ 



J2 



d'9 up{92)ug{9i - 92)un-p-q{9 - 9i) 



(22) 

which is an integral equation for the set {un}nez of 
Fourier coefficients. The depinning exponent ()11[) in the 
Fourier space becomes 



B{T >0,h) = 




UniSl)Umi92)U-n^rr.{-9i - 92) - ^ V / d^^" 
Un{9l)Um{92)ui{9z)U-n-ra-l{-9l~92~9:i) Q (23) 

The numerical algorithm is analogous to the one used 
in the zero-temperature case, but taking into account 
the reescaling of {-Upjpgz by a factor -\/27rA„/3wo(/i) and 
{Up}p^z by a factor V27ra„ in the last step of the calcu- 
lations. 



d^e 
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FIG. 2: (Color online) Temperature dependence of the de- 
pinning rate; normalized action - — ^^^=_B(T) versus S7 at 



|pg| V kA 



different values of the parameter h. 

Fig. [5] shows the normalized action B{T) — 
^ -B{T) as a function of 51 at different values of 
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the parameter h. In the simulations we have taken the 
standard value a^Lc — 0.008 for bulk Permalloj^. 



IV. CROSSOVER TEMPERATURE 

The crossover temperature determines the transition 
from thermal to quantum tunneling relaxation regimes. 
It can be computed by means of theory of phase 
transitions^: above Tc, the instanton solution minimiz- 
ing Eq. PT|) is a f-independent function u{f, z, h) = 
UQ{z,h), whereas just below the instanton solution 
can be split into the sum of uq and a small perturbation 
depending on f, 
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27r 

u(f , z, h) — uq{z^ h) + ui(z, h) cos ( — f 



The depinning exponent (jlip is proportional to 



/ dz 



(24) 



(25) 



where $ is the spatial action density. Introducing the 
expansion ([M)) into Eq. ([TT]) we obtain the following 
expansion 



n 



(u'i)^+Auf + 0(4) 



with 



(26) 



FIG. 3: (Color online) Field dependence of the dimensionless 

47rfc_g I PG I 



crossover temperature Oc = 
the potential V{u, h — 0.1). 



'-Tc. (Inset) Sketch of 



with f{h)~4: — 2ul{h). Solving the quadratic equation 
for Tc given by Eq. ([SS)) we obtain the crossover temper- 
ature 



Tcih) = 



47rfcB I pG I 



l8 + 3f{h)-12uo{h)./f(h)+ 2 11 

Pg \pg\ 

(31) 

Figure [3] shows the dependence of the dimensionless 
ossover temperature 
ized magnetic field h. 



crossover temperature 9c = b\pGj_j^^ general 

Hk 



A = Sv>,»4 + ^^ (27) 

If A > the only pair (ui.,Ui) minimizing $ is ui = 0. 
The crossover temperature is then defined by the equa- 
tion min^gK A = 0, that is 



^ min V"{uo. h) + — + ^j^ 







(28) 



The equation of motion for a r-independent instanton is 

Uq - (2 - 3ul{h))uo + 3uo{h)ul + = (29) 

with boundary conditions: mq — > at |z| — > oo and 
Mo(0, h) = -2uo{h) + a/4 - 2u1{h) = w{h), which is the 
width of the potential. Consequently, 



mmV''(uo(z,h)^h) = min V"(uQ,h) 

zm uoelo,w{h)] 



mm 

uoe[o,w{h)] 



[{2 - 3ul{h)) - 6uoih)uo - Ml'j 



10 + 3ul{h) + 6uo{h)./J(li) 



(30) 



V. DISCUSSION AND PARAMETERS FITTING 

For a given value of the generalized field h, in Fig. [2] 
we clearly distinguish two regimes in the dependence of 
the normalized action on fi: above ildh) the normalized 
action tends to a constant value, whereas below it the 
normalized action is linear with fl. Notice that the tran- 
sition from the linear to the constant regime is smooth 
(that is, of second-order type). Above Tc the depinning 
rate becomes 



l{u'of + Viuo,h) 



n 

(32) 



with Uq being the r-independent instanton. By means of 
Eq. ([29| this expression can be rewritten aa^ 



B{T>Tc,h) 



2hl3 



aih) 



2V2 / duoVV{uo,h) 



n 

(33) 



and, consequently, the slope of the normalized action 
B{n) is equal to 2\/2 J^^''^ duoy/V{un, h), which can be 
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evaluated analytically. At all values of the generalized 
field h, the numerical slope calculated from Fig. [5] coin- 
cides with the analytical one within the numerical error 
of our simulations. This is indicative of the robustness of 
our algorithm. 

Quantum effects reported in Ref. HO can be understood 
as being plausibly due to the depinning from line defects 
present in the disk. The size of the defects needs to ex- 
ceed the nucleation length in order to pin the VC, but 
not to be as long as the thickness of the disk. Pinning 
of extended parts of the VC line by line defects would 
be justified by the fact that linear defects provide the 
strongest pinning so that the VC line, or at least some 
segments of it, would naturally fall into such traps. Con- 
sequently, we can test out our model on the experimental 
results obtained in Ref. ,20, . The crossover temperature is 
relevant to the roughness of the fine-scale potential land- 
scape due to linear defects at the bottom of the potential 
well created by the external and dipolar fields. Above 
Tc vortices diffuse in this potential by thermal activa- 
tion, whereas below Tc they diffuse by quantum tunnel- 
ing. This must determine the temperature dependence 
(independence) of the magnetic viscosity. Tc is, there- 
fore, the measure of the fine-scale barriers due to linear 
defects. It can be measured experimentally and help to 
extract the width of the pinning potential. 

Now we proceed to obtain estimates of the model pa- 
rameters (k, /3) by fitting our model to experimental data: 
Figure U shows new magnetic relaxation measurements 
of permalloy disks in the vortex state from the remnant 
state to equilibrium (zero magnetization). The radius of 
these disks is R = 0.75 fim and their thickness is L = 95 
nm (subfig. |4^) and L = 60 nm (subfig. Hb). A con- 
cise description of the experimental set-up and sample 
preparation can be found in Ref. HO- Notice that for 
both samples the magnetization depends logarithmically 
on time during the relaxation process. 

Magnetic viscosity of these relaxation measurements is 
computed by means of the formulai^l: 



S{T) 



I dM 



(34) 



where Mq is the initial magnetization point. That is, the 
viscosity at zero field is obtained computing the slopes 
of the normalized magnetization curves. Fig [5] shows the 
magnetic viscosity as a function of temperature for both 
samples. Below Tc — 6 K, magnetic viscosity reaches a 
plateau with non-zero value. Above Tc, magnetic viscos- 
ity increases up to a certain temperature, from which it 
decreases again. The existence of the plateau is the ev- 
idence of underbarrier quantum tunneling phenomena. 
The increase of viscosity with temperature above the 
crossover temperature is due to thermal activation over 
the pinning barriers. Finally, the drop of the magnetic 
viscosity is in agreement with the loss of magnetic irre- 
versibility in our systems^S. On the other hand, the fact 
that the crossover temperature Tc is independent of the 
thickness of the disks upholds our hypothesis that just a 
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FIG. 4: (Color online) Relaxation measurements of magnetic 
vortices from the remnant state to equilibrium for samples 
{L,R) = (95,750) nm (subfig (a)) and {L,R) = (60,750) nm 
(subfig (b)). 
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small portion of the VC line takes part in the tunneling 
process via an elastic deformation. 

Notice that the depinning rate should not exceed 30 — 
40 in order for the tunneling to occur on a reasonable 
time scale. The estimates of the parameters (k, /3) are 
obtained fitting Eq. ([31]) and Eq. (|16p to the values 
Tc ~ 6 K, respectively B{T = 0, /i = 0) - 30 at zero field. 
Considering the experimental values A = 1.3 • 10~^^ J/m 
and Ms — 7.5 • 10^ A/m for permalloy, we obtain 

K - 5.9- 10^ J/m^ /3 - 6.9 • 10^^ J/m^ (35) 

from which we can determine the width of the quartic 
potential, w — \/2k/ f3 ^ 0.13 nm. This value is compat- 
ible with the width of the potential provided by a linear 
dislocation 

In conclusion, we have studied quantum escape from 
a line defect of the VC line in a disk made of a soft 
ferromagnetic material. In the case of permalloy disks. 



experimental results let us conclude that the depinning 
process occurs in steps about 0.13 nm, which corresponds 
to the width of the energy potential. 
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